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Three-dimensional nonlinear finite-element heat transfer and structural 
analyses were performed for the first stage high-pressure fuel turbopump blade 
of the space shuttle main engine (SSME). Directionally solidified (DS) 

MAR-M 246 material properties were considered for the analyses. Analytical 
conditions were based on a typical test stand engine cycle. Blade temperature 
and stress-strain histories were calculated using MARC finite-element computer 
code (ref. 1). This study was undertaken to assess the structural response of 
an SSME turbine blade and to gain greater understanding of blade damage mecha- 
nisms, convective cooling effects, and the thermal-mechanical effects. 

Hot-gas path components for reusable space propulsion systems operate 
under extreme gas pressure and temperature. These operating conditions subject 
the high pressure stage turbine nozzles and blades to severe thermal transients 
that can result In large Inelastic strains and rapid crack Initiation. 

Advances In casting techniques have allowed the development of directionally 
solidified and single crystal alloys for high temperature components In space 
propulsion vehicles. Mechanical anisotropy exhibited by these alloys have to 
be taken Into account In the analytical studies. To Improve the durability and 
accuracy of turbine blades and other hot section components, an accurate knowl- 
edge of the temperature and stress-strain histories at the critical location 
for crack Initiation Is required (fig. 1). 

Experimental measurements of gas temperature profiles In the SSME turbo- 
pumps are difficult to obtain because of the high gas temperature and severe 
thermal transients. Turbine blade temperatures are primarily a function of hot 
gas flow and cooling. The temperature field Is determined by the heat transfer 
from the hot gas to the blade. This heat transfer and Its variations are 
determined by knowledge of the gas film coefficients. Also, the time- 
temperature history profile at the start transient, steady state, and cutoff 
Is obtained through a combination of analytical and experimental results which 
Is due to the complex flow phenomena through an accelerating turbine. 

Temperature-dependent properties for the MAR-M 246 + Hf alloy were mainly 
provided by Rockwell International Corporation. These elastic properties are 
summarized In table I. Mean thermal coefficient of expansion data were con 
verted to Instantaneous values for MARC Input. Longitudinal stress-strain 
properties, summarized In table II, were used for elastic-plastic region. The 
mission used for this analysis Is shown In figures 2 and 3 In terms of Inlet 
temperature, gas pressure, and revolutions per minute (RPM). This cycle Is 
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applicable to a factory test of the engine; It Is also reasonably representa- 
tive of a flight mission except for the foreshortened steady-state operating 
time. The major factor Inducing fatigue cracking Is the transient thermal 
stresses caused by the sharp Ignition and shutoff transients. The finite- 
element model of the blade In figure 4 was constructed of a three-dimensional 
eight-node Isoparametric brick element. The model consisted of 1-25 elements 
with 1575 nodes and 4660 unsuppressed degrees of freedom. 

Heat transfer coefficients at the blade airfoil (figs. 5(a) and (b)) were 
predicted by running a boundary layer analysis using a modified version of the 
STAN5 boundary layer code (refs. 2 and 3). Thermal environment experienced by 
the platform and shank were obtained from reference 4. Details regarding pre- 
diction of heat transfer coefficients at the stagnation region for the airfoil 
are available In reference 5. Transient results were obtained by scaling 
steady-state heat transfer coefficients based on transient flow and tempera- 
ture. Predicted high temperature locations were evaluated for the two tempera- 
ture spikes shown In figure 2. 

The thermal response predicted from the finite-element analysis showed 
that the leading and trailing edges of the airfoil base are the hottest loca- 
tions. Temperature distributions showed a cordwlse variation at the first 
Ignition spike and a spanwlse variation thereafter Into the cycle. A uniform 
temperature distribution was dominant at most of the airfoil surface during 
cruise except near the base at the platform junction where a mixture of cold 
and hot gas Is present. The coldest spot was always at the blade root because 
of cooler boundary conditions. Blade temperature Is shown In figures 6 
through 8. 

Elastic-plastic analyses have been conducted for the HPFTB blade with MARC 
code. Plastic strain calculations were based on Incremental plasticity theory 
using Von Mlses Yield criterion, the normality rule and a kinematic hardening 
model. The material elastic-plastic behavior was specified by the yield 
strengths and work hardening properties In the longitudinal direction; trans- 
verse properties were not available. Creep analyses were not performed at the 
present time because of Inadequate knowledge of the creep characteristics for 
anisotropic blade material. 

Incremental loading Included centrifugal and gas pressure loads and metal 
temperature distributions as calculated from the heat transfer analysis. The 
same Increments were used for the heat transfer and the elastic-plastic struc- 
tural analysis. Approximately two million words of core storage on the CRAY- 
XMP computer were needed to run the problem. Analysis required about 4 hours 
of central processor unit (CPU) time on the CRAY system. The directionality 
of the elastic material properties causes anisotropic constraints. Lekhnltskll 
(ref. 6) has derived the generalized elastic strain equations for an aniso- 
tropic body with a transverse plane of Isotropy. This anisotropic stress- 
strain law was Incorporated In MARC user subroutine HOOKLW. 

The colder airfoil base temperatures Induce tensile thermal stresses at 
the critical leading edge location that are additive to the centrifugal 
stresses. A noticeable discrepancy Is seen In the compressive strain region 
where the reduction In elastic analysis points resulted In failure to capture 
some of the cycle fluctuation due to transient thermal effects during the rapid 
engine cooldown. The CPU time for the elastic finite-element analyses amounted 
to 5 percent of that required for one cycle of the nonlinear finite-element 
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analysis. The critical location for crack Initiation Is at the leading edge 
near the base. Results of both elastic and Inelastic structural analyses Indi- 
cated the region of the finite-element model with the largest total strain 
range was coincident with the observed crack Initiation site. 

Radial stress distributions as calculated at the gausslan Integration 
points closest to the suction and pressure surfaces are presented In figures 9 
through 11. The stresses primarily reflect the centrifugal and thermal load- 
ings. It Is apparent that the analysis would not show a low-cycle fatigue 
problem If transient thermal effects were not considered. 
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TABLE I. - OS MAR-M 2*6 PHYSICAL PROPERTIES 


Temperature, 

•c 

Modulus of elasticity, 
GPa 

Means coefficient of 
thermal expansion, 
percent/*C 

Longitudinal 

Transverse 

21 1 

jHRHH 
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204 


175 

.00130 

316 


173 


427 

119 

166 

.00141 

538 

114 

162 

.00148 

649 

109 

156 

.00149 

760 

103 

149 

.00156 

671 

97 

142 

.00160 


TABLE II. - DS MAR-M 246 STRESS- 
STRAIN PROPERTIES 
(LONGITUOINAL) 


plastic strain, 
percent 


0.1 

.2 

.4 

.6 

.6 

1.0 


Stress, MPa 

21 # C 


B9BI 

600 


875 

B30 


930 

850 


965 

855 

930 

970 

865 

945 

975 

870 

960 

980 
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Figure 1. - SSME high-pressure fuel 
turbopump first stage turbine blade. 
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(a) Startup transient. 


(b) Cutoff transient 


Figure 2. - Turbine Inlet gas temperature for HPFTP. 
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(a) At airfoil root. 
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(b) At airfoil tip. 

Figure 5. - Distribution of heat transfer coefficients. 
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